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The unexplained tunnehng conductances of correlated mesoscopic Kondo systems are understood 
by the coherent transport of the entangled and nonentangled singlets. Spins of the entangled singlet 
flow unidirectionally in a sequential up-and-down manner. This dynamics does not follow linear 
response theory. The side peaks at a finite bias are formed by resonant tunneling of the nonentangled 
singlet through a coherent tunneling level formed by two electron reservoirs within a coherent region. 
The theoretical line shapes remarkably fit the experimental data of a quantum point contact and a 
magnetized atom adsorbed on an insulating layer covering metallic substrate. 



The nonlinear line shapes of the tunneling conductance 
observed for mesoscopic Kondo systems are waiting for 
relevant theoretical explanations from the microscopic 
point of view. The two-reservoir Anderson impurity 
model at steady-state nonequilibrium is considered as a 
proper microscopic model describing a mesoscopic Kondo 
system. However, previous theoretical studies using 
the noncrossing approximation p]], the Keldysh formal- 
ism quantum Monte Carlo calculations [3], and the 
extended numerical renormalization group method [4] do 
not reproduce various nonlinear line shapes of mesoscopic 
Kondo systems. The difficulty originates from the combi- 
nation of strong correlation and steady-state nonequilib- 
rium. To solve this problem, we need to understand the 
dynamics under bias more clearly. It is obvious that the 
two-reservoir mesoscopic Kondo system at equilibrium 
has an entangled Kondo singlet represented by a wave 
function |^)en = p\i't)LS + q\U)LS + u\it)RS + v\U)rs^ 
where LS and RS denote left singlet and right singlet, 
respectively and the coefficients are complex numbers, as 
shown in Fig. 1 (a). The equilibrium dynamics comprises 
multiple processes of exchange, partner change, and sin- 
glet hopping in a mixed manner among the four states 
given above. This complicated dynamics may be studied 
using the numerical renormalization group method for 
the low-energy regime [5]. Under bias, however, the dy- 
namics becomes considerably simpler because the spins 
involved in the entanglement flow unidirectionally in an 
up-and-down sequence, as shown in Fig. 1 (b), until the 
entanglement is retained. The processes of singlet hop- 
ping and partner changing are used in the spin flow, and 
the coherent spins are provided from the Kondo cloud. 
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FIG. 1: (Color online) (a) An entangled Kondo singlet at 
V = 0. The dashed circle and the letter "m" denote the 
Kondo cloud and the mediating Kondo atom, respectively, 
(b) An entangled Kondo singlet at V ^ 0. Singlet hopping 
and partner changing leads to unidirectional spin flow. The 
horizontal arrow denotes the direction of the electric field. 



The unidirectional fiow does not allow a linear response 
regime. Backward motion of electrons may occur in the 
incoherent dynamics, and this leads to double occupancy 
at the mediating atom. 

In this study, we clarify the spin dynamics forming 
the zero-bias peak and the side peaks and compare the 
theoretical results with the experimental data obtained 
for a quantum point contact [6] and an adsorbed mag- 
netized atom on an insulating layer covering a metallic 
substrate (7], for example. Fitting the entire range of the 
line shapes of these systems is given for the first time. 

In Fig. 2, we depict the low-energy tunneling schemes 
in the two-reservoir Anderson impurity model under bias 
along with the corresponding line shape of the tunnel- 
ing conductance. Figure 2 (a) describes the transport 
near zero bias, i.e., eV < ksTx^ where e, V, /cb, and 
Tk denote the electron charge, source-drain bias, Boltz- 
mann's constant, and Kondo temperature for the entan- 
gled Kondo singlet, respectively. The unidirectional fiow 
described in Fig. 1 (b) establishes the zero-bias peak in 
the tunneling conductance. The entanglement is com- 
pletely broken when eV > 2kBTK^ in which a nonen- 
tangled Kondo singlet, \^)ne = pIWls + qI^Dls, per- 
forms resonant tunneling and yields the side peak when 
it reaches the coherent transport channel shown by the 
dashed line in Fig. 2 (b). This transport also prohibits 
backward motion. Therefore, prohibiting backward mo- 
tion of coherent spins is a generic feature of transport 
in a mesoscopic Kondo system under bias. This prop- 
erty significantly simplifies the dynamics at steady-state 
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FIG. 2: (Color online) (a) Transport of an entangled Kondo 
singlet at low bias, (b) Resonant tunneling of a nonentangled 
(N. E.) Kondo singlet through the coherent tunneling level 
(dashed Une). Sm denotes the energy level of the mediating 
spin that forms a singlet. 
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nonequilibrium. Another crucial feature is the existence 
of two coherent transport channels, shown in Fig. 2 (b), 
which is attributable to two reservoirs within the coher- 
ent region. Specific proof is given in a previous study Q 
and the basis vectors given later in the text clarify their 
existence. Previous studies [H-Q have obtained a similar 
result showing Kondo peak splitting with bias. This phe- 
nomenon may occur when the two reservoirs are out of 
coherence. However, each of the aforementioned meso- 
scopic systems should be considered as a complete coher- 
ent system. 

Now, we validate the tunneling mechanisms given in 
Fig. 2 by obtaining the tunneling conductance that fits 
the experimental result. The tunneling current of a meso- 
scopic system with an interacting site between two non- 
interacting reservoirs is given by [9j 



where fL{R){^) is the Fermi distribution function of the 

left (or right) reservoir, f{uj) = T^{uj)T^{uj)/[T^{uj) + 

T^{uj)i where r^(^)(a;) = 27rZk l^fcrn""^ - ^fe) ^^at 
involves the reservoir density of states, and pf^{oj) is the 
local density of states (LDOS) at the mediating atom. 
The superscript ss means steady-state nonequilibrium. 
This simplified form is derived from the well-known Meir- 
Wingreen current formula [10, 11] by using the condition 
of the proportionate coupling function, T^{uj) ex r^{u). 
This proportionate relation can be applicable to the elas- 
tic tunneling shown in Fig. 2. We employ a constant 
r^'^(a;), indicating a fiat density of states of metallic 



reservoirs. 



Since the electrons in the singlet do not collide with 
a quasiparticle such as phonon until the bias excites 
it, p^{uj) in this study must be bias independent, i.e., 
dpf^{uj)/dV = 0, which means one can write the tunnel- 
ing conductance at zero temperature as 

dl/dV = {e/h)f{oj)p^^{uj)\n^=,v 

The bias independence of p^{uj) will be retained unless 
inelastic tunneling is caused by the scattering with quasi- 
particles. In other studies on Kondo-involved mesoscopic 
systems [12-14], dpf^{uj)/dV has also been neglected. 

We obtain pf^^{uj)^ which is given by p^^{uj) = 
— (l/7r)Im[Gj^^^(co')], by calculating the on-site retarded 
Green's function, iG^^^^{uj) = (c^tK^ + zL)"^|c^^), 
where z = — icj + O^, Cm^ is the fermion operator that an- 
nihilates an up-spin electron at the mediating atom and 
L is the Liouville operator defined by LA = HA — AH^ 
in which H is the Hamiltonian and A is an operator. 
To obtain G^^^(co') using the resolvent form, one needs 
a complete set of basis vectors spanning the Liouville 
space. We have obtained a complete set of ort honor mal 
basis vectors [H, [l5| describing Cm^{t) that is driven by 
the Hamiltonian of the two-reservoir Anderson impurity 



model. 
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^Uflmt^mi^y^C, (1) 



where Ti 



L,R _ 



'^k,a,u^L,R(^km^tia^ka + "^/cm^/ccr^^cr) , and (J, 6/^, e^n, 

Vkm^ U and p indicate the electron spin, kinetic en- 
ergy, energy level of the mediating atom, hybridization 
strength, on-site Coulomb repulsion, and chemical poten- 
tial, respectively. 

We divide the complete set of basis vectors at equilib- 
rium [8] into four groups: 

H: {c^j , nm^c^j , J^^c^-f ? Jmf^kt I ^ = 1' 2, • • • , oo}, 
HI: {{L^jt'ty^^, (^cJ™f )cmt I n = 1, . . . , oo}, 



±R\^R 



n, 



1,2, 



, oo 



where Lc denotes the Liouville operator using Hc^ 

Jmi = Eki'^kmcl^^Cki + Vfc*m4iCmi), and j-^ = 

^E/c(^fcmc]^;C/e; - Vj^^cl^Cmi)- Groups HI and IV rep- 
resent multiple trips of a down-spin electron, and they 
do not play a role in describing the unidirectional mo- 
tion of the Kondo singlet. As a result, the degrees of 
freedom of the system are considerably reduced when a 
bias is applied. The degrees of freedom are further re- 
duced by neglecting the basis vectors rimiCmt group I 
and j^^Ck^ in group II. The former describes all higher 
orders of double occupancy, i.e., from U to If^^ because 
^mi ~ ^rni- Therefore, the basis vector rimiCm^ must 
be neglected in this study. In contrast, the latter basis 
vectors are neglected because of their minor contribution 
to self-energy compared with UmiCk^ in group II, whose 
members play the role of constructing self-energy. Now, 
we have a working Liouville space [8] that is spanned 

by {^imt^rnt, ^J^t^mt. C^t^ ^J^f^mt. ^J^f^rnt} and 

{c^^, Srimic^^, SrimiC^^, c^^\k = 1,2, •••,00}. We use 
5 indicating 3 A = A — (A), where the angular brackets 
denote the expectation value, to achieve orthogonality 
among the basis vectors. For convenience, we omit the 
normalization factors '^)^)"^^^ and {{Srimi)'^)^^'^ in 

the denominators of the corresponding basis vectors. 

We construct the matrix Moo x 00 = zI-\-iL in terms of 
the basis vectors spanning the working Liouville space. 
Matrix reduction from M 00x00 to Msxs is possible for 
noninter acting reservoirs [16|. Thus, p^^{uj) is given by 



using the relation p^.(co') = (l/7r)Re[(M5 



5x5; 
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Further, 
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where 



(2) 



= uj - Cm-U (rimi) and (n^;) denotes the 
average number of down-spin electrons occupying the me- 
diating atom. All the matrix elements, except U^±^ ^ 
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have additional self-energy terms iSrnn = Pmn[^^o i^) + 
where ^o^^\uj) = -iT^^^^ /2 for a flat wide 
band. We use A = (F^ + r^)/4 as an energy unit. The 
coefficients /3mn are discussed in the following text. 

The matrix in Eq. ([2]) consists of two 3x3 blocks 
that share the central element representing the mediating 
atom and two 2x2 blocks at the corners. The three off- 
diagonal elements of the 3x3 block represent the degrees 
of singlet coupling J^^^^^y e.g., 



1l 



and the incoherent double occupancy parameters 

■±L,R\ 



where (Z)±^f ) = (K;, - 2n^,))/(i±f ^). 
Hence, the real part Re[U^±^] represents the probabil- 
ity of double occupancy by or coming from the 
left or right reservoir. In contrast, the elements of the 
2x2 corner blocks, i.e.. 



^km ^fct ) ^^mi ' ^mf ] ) 



and 7 ^ 



7 . , where 



{EM 



km^kt 



V* r 



1/2 



describe the transition between two reservoirs. In Fig. 3, 
we present graphical illustrations of the third order hy- 
bridization processes embedded in 7 under bias. Singlet 
partner change and singlet hopping will occur to perform 
the process of Fig. 2 (a). In contrast, only singlet hopping 
is needed for Fig. 2 (b). The basis vectors in groups III 
and IV cause incoherent or backward motion that must 
be excluded in describing the unidirectional coherent mo- 
tion at steady-state nonequilibrium. Therefore, neglect- 
ing groups III and IV is legitimate. The unidirectional 
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FIG. 3: (Color online) Spin dynamics under bias in 7. The 
numbers denote the sequence of coherent motion. A down 
spin moves into the mediating atom (blue circle) and performs 
an exchange [green 2 in (a)] in 7^^ and a singlet hopping 
[red 2 in (b)] in J^ji j- The second part of (b) vanishes because 
of the reverse motion. 



movement of the Kondo singlet discussed in Fig. 2 guar- 
antees 7^^ = J. in Fig. 3 (b). This equality is considered 
as a condition of steady-state nonequilibrium. 

The 5x5 matrix of Eq. ([2]) gives three coherent and 
two incoherent poles in the LDOS. One of three coher- 
ent poles is located at the Fermi level and the other 
two are at ±hujrti^ which are the levels represented by 
the dashed lines in Fig. 2 (b). The subscript rti indi- 
cates the resonant tunneling level. Hence, the tunneling 
current rapidly increases when the bias voltage reaches 
:^hujrt£- A simple atomic limit analysis using the same 
Rep^''] gives hcoru = ±[hl, + jljm^/^ + 0{U-^) 
and the spectral weight of the zero-bias peak as Zjr = 



lLllJ2Re[U^^]H^l^ + 7L) + OiU-^). The latter 



,,,ylJ2Re[U^^niL ■ -KH. 
expression shows that the zero-bias peak is suppressed 
when 7^^ and 7^^ are imbalanced. 

We show in the following that Re[U^^^] > Re[U^^^] 
is introduced for symmetric reservoirs, which means that 
operator induces more double occupancy than 
For asymmetric reservoirs such as substrate (L) and tip 
{R) in scanning tunneling spectroscopy (STS), we use 
Re[U^+] > Re[[/^_'^] > Re[U^+] to reflect different prop- 
erties of reservoirs. In contrast, we adopt Re[U^_] = 
Re[U^] because of the steady-state property of the cur- 
rent operator 

Before obtaining the dl/dV line shapes, we first deter- 
mine the coefficients Prnn- ^^[^25]^ for example, is given 
by Re[/?25] = - 2n™t))0-+f (1 - 2n^t)) + (1 " 



2(n™,))^0lt)0™f)}/4{(<5n„i)^)V^((5j+t)^)v/((^i™f)^)- 
The operators in the first term describe the self-energy 
dynamics that avoids double occupancy Q. We as- 
sume the same contribution of self-energy dynamics 
to all Re[/3mn] and the same relative fluctuations. 
Therefore, the differences among Re[/3mn] are at- 
tributable to the different signs for the current, i.e., 
iJ^i) = -0™f) < 0. From the property > 0, 

symmetric Re[/3^n] have the following mutual relations: 
Re[/3i2] = Re[/3i4] = Re[l3,,] < Re[/3n] = Re[/322] = 
Re[/344] = Re[/355] = Re[/324] = Re[p25] = MP^^] and 
Re[/333] = 1. In this study, we choose Re[/3ii] = 0.255 
and Re[/3i4] = 0.245 based on the standard value 
Re[/3mn] = 0.25 that is obtained at the atomic limit [H. 
The same relative fluctuations give vanishing Im/^^n 
because they are given by the difference in the relative 
fluctuations. 

The experimental dl/dV line shapes under considera- 
tion are those of a quantum point contact with the closest 
side peaks given in Fig. 1 (b) of Ref. [6] and the STS for a 
Co atom placed on a CU2N layer on a Cu (100) substrate 
of Ref. [7J. For the former, we employ the scenario of 
spontaneous formation of a localized spin at the bound 
state [H, [i^. Therefore, the Hamiltonian of Eq. ([!]) is 
applicable to both cases. The gate voltage dependence 
in the former system will be described in a future study. 
The theoretical dl/dV line shapes given in Fig. 4 are ob- 
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FIG. 4: (Color online) Comparison of the theoretical line 
shape (red) using A ^ 1.5 meV (a) and A ^ 5 meV (b) 
with the experimental data (black) of the closest side peaks 
reported in Ref. [6] and the STS line shape of Ref. ^7,] , respec- 
tively. We choose F/A = 0.82 in (a) and an arbitrar y u nit for 
the theoretical dl/dV in (b). Inset is STS setup f?,\l§\. 



Table I: Matrix elements for Fig. 4 
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^j,LR 


Re[/^'^ 
J 




ReUf+ 


3^ 


(a) 


0.38 


0.38 


0.5 


0.98 


1.1 


1.1 





(b) 


0.86 


0.76 


0.43 


2.8 


7.0 


1.62 






tained by using the matrix elements given in Table I. We 
set the energy unit to A ^ 1.5 and 5 meV for Figs. 4 
(a) and 4 (b), respectively. The fittings are remarkably 
good. 

As listed in Table I, we adopt a symmetric Kondo cou- 
pling (7,, = 7,,) and Re[C//;] = Re[C7f^] > Re[U^l''] 
in Fig. 4 (a) and an asymmetric Kondo coupling (7^^ ^ 
7,,) and MUfA > Ke[/7/_'^] > Re[/7f+] in Fig. 4 (b). 
The zero-bias peak in Fig. 4 (a) is slightly suppressed by 
different contributions to double occupancy by operators 
j+ and j~ and the position of the side peak matches huru 



pretty well. The deviation outside the side peaks indi- 
cates that the range of bias independence of the LDOS 
covers the two side peaks. It is noteworthy that the line 
shape of Fig. 4 (b) shows a dip at zero bias when the in- 
sulating layer is removed [17]. This implies that a strong 
Kondo couphng is established by inserting the insulat- 
ing layer. Choi et al. [l8^, who studied the same system, 
observed that the meaningful structure of the line shape 
disappears when a Co atom is placed on top of a N atom. 
This indicates that the N atoms surrounding a Cu atom 
in a CU2N layer play the role of barrier that suppresses 
fluctuations and enhance Re[U^+] and the axial Kondo 
coupling connecting tip, Co atom, and Cu substrate. The 
large values of , and given in Table I verify this 
fact. Choi et al. [18] also show that the shoulder in Fig. 4 
(b) is a variation of a coherent side peak. 

In conclusion, our theoretical study clarifies that there 
are two different transport channels (Fig. 2). One uses 
the entangled Kondo singlet that connects the Kondo 
clouds in both reservoirs. Transport by the entangled 
Kondo singlet forms the zero-bias peak. The other 
uses resonant tunneling of a nonentangled Kondo sin- 
glet through the coherent tunneling level. This tunnel- 
ing mechanism forms the side peak. Comparisons of the 
theoretical dl/dV line shapes with those of the experi- 
mental ones (Fig. 4) clearly demonstrate the existence of 
the two transport channels. 
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